NON-PERTURBATIVE ASPECTS OF MULTIPARTICLE PRODUCTION 



in 

0\ V.A.RUBAKOV 

Institute for Nuclear Research of the Russian Academy of Sciences, 60th October Prospect la, 
> . Moscow, 111312, Russia 

O ' E-mail: rubakov@ms2.inr.ac.ru 

Iz; : \hev-vh/mi2m 



> 

m 
(N 

ABSTRACT 

in 

On \ Processes with multiparticle final states in weakly coupled theories, both with 

and without instantons, cannot be studied perturbatively at most interesting 
^ i' energies and multiplicities. Semiclassical approaches to the calculation of their 

total cross sections are reviewed. 

^ ■ 1. Introduction 

X, 

^ , In this paper I discuss interacting bosonic theories (fermions are beheved to play minor role 

in the processes to be considered) whose coupling constant, generically denoted by a, is small. 
At cm. energies of order m/a, where m is the typical mass of elementary quanta, the number of 
particles kinematically allowed to be produced is of order l/a. It is this multiplicity that will be 
of interest in this paper. The non-perturbative character of the the processes where the number 
of outgoing particles is so large may be seen already from combinatorics of the corresponding 
graphs. On the other hand, energies of interst are not exponentially high, « In ^ -C 1, therefore, 
such aspects as running of the coupling constant or triviality are irrelevant. 

There are two classes of multi-particle processes that are discussed most widely. Historically, 
considerable interest has been attracted by instanton-like processes 0i which, in the electroweak 
theory, violate baryon and lepton numbers i. It has been realized Bi that multiparticle pro- 
cesses without instantons show remarkably similar behavior. From purely theoretical point of 
view, both classes of processes are of interest because leading order calculations indicate rapid 
increase of the cross sections at high energies and multiplicities, the "corrections" are more 
than important at energies and multiplicities of order 1/a and, ultimately, novel techniques of 
semiclassical type should be developed for the complete treatment of the cross sections. On 
the phenomenological side, the processes under study include multiple W- and Higgs boson 
production (with or without violation of B and L) in the electroweak theory at energies in 



multi-TeV range and also possibly multiple production of semi-hard gluons in QCD; clearly, 
these phenomena, if occur at reasonably high rates, are of general interest. 



2. Perturbative Results 

Let me briefly summarize perturbative results on the multi-particle processes; in the case of 
instanton-like processes, "perturbation theory" means perturbation theory about an instanton. 
I will not discuss most of perturbative techniques, even though they are of interest by themselves 
and provide some important hints for non-perturbative analysis; for reviews see refs. iSi. 

2.1. Instanton- Like Processes 

As a prototype model for instanton-like processes at high energies, consider an SU (2) gauge 
theory with one Higgs doublet (minimal standard model in the limit sin^ 6'^^- = 0). In the 
leading order of perturbation theory about an instanton, the cross section of the oroduction of 
a given number of W- and Higgs bosons rapidly increases with cm. energy E aB. The total 
instanton-like cross section in the leading order shows the exponential behavior iSEUll, 



leadinq order 
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where Eq = y/Enmw / dw is of order 15 TeV in the electroweak theory. The number of IV-bosons 
produced is of order 

An fE\^l^ 

n 



aw \Eq, 

The cross section (|l]), being exponentially small at low energy, exponentially increases with 
energy and hits the unitarity bound at E Eq; at these energies the multiplicity of outgoing 
particles becomes of order n ~ 1/aw 

The expression ([^) should not determine the true cross section at all energies, as it violates 
unitarity. Indeed, corrections to the leading order result become large aX E Eq. Generally, 
the total cross section has the following functional form O 
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and perturbation theory about the instanton provides the expansion of F{E/Eq) in 

powers of (E/Eq)^/^, 

\eJ 8 \EoJ 16 \EoJ 

This series blows up at ~ i?o, and new non-perturbative techniques are needed to find the 
exponent of the cross section, F, in the interesting energy range. 

A few remarks are in order. First, the exponent F{E/Eq) in Eq. (|^) receives contributions 
both from tree graphs (about an instanton) and from loop graphs; the latter involve hard 
incoming particles tlS with energies of order E. Second, Eq. (|^) looks semiclassically, and 
one naturally expects that there exisis semiclassical technique for calculating the exponent. 
This semiclassical technique should incorporate loop effects about an instanton (!). Third, the 
behavior of F[E/Eq) (whether it becomes zero at some energy or remains always negative) 



determines whether the total instanton cross section becomes unsuppressed at some energy 
or is exponentially small at all energies. Finally, the general form of Eq. (|^) is valid in all 
weakly coupled models with instantons, though the perturbative series for F{E /Eq) is model- 
dependent. 



2.2. Processes without Instantons 



As a prototype model with multi-particle production and trivial vacuum consider massive 
theory without spontaneous breaking of the symmetry (j) —cp. The direct analysis of tree 
diagrams with n outgoing particles is possible at exact threshold, i.e., when all final particles 
are at rest. The exact result for the tree amplitude at threshold is 

A^-SE = nm) = n\ (^^j ' (3) 

The factorial behavior of the amplitudes indicates that the tree cross sections also increase with 
n (if amplitudes do not drop too rapidly with energy) , 

~ ^l4Z'„r X {phase space) ~ n!A"e" (4) 



where 



[E — nm)/n 



is the average kinetic energy of outgoing particles. Indeed, the direct analysis of the tree 
diagrams ils has lead to a lower bound for the tree cross section which grows factorially with 
n in a way similar to Eq. @. 

At n ~ 1/A the tree cross section grows with n and exceeds the unitarity limit at 
large enough n; the corresponding total cm. energy is of order m/A. This means that loop 
"corrections" must be large at these multiplicities. In fact, the one loop contribution @ to the 
tree amplitude at thershold becomes large even before the multiplicity becomes of order 1/A, 

+ K:'r = + sxn^) (5) 

where (complex) numerical coefficient B is of order one and has been calculated in refs. 
Further perturbative analysis"has shown HI that both the loop corrections and energy correc- 
tions exponentiate, which provides a strong argument showing that the true cross section has 
exponential behavior, 

ai^n{.E) oc exp [nF{Xn, e)] (6) 

in the regime of small A, large n and An, e = fixed. The loop expansion for the exponent F is 
the expansion in (An); only a few terms in the expansion of F at small An and e are known, 

f(An,.) = l„- + - + -l„---e + BA„ + ... (7) 

In complete analogy to the instanton case, the cross section (^) shows the semiclassical behavior, 
and the perturbation theory for the exponent F blows up in the most interesting region An ~ 1 
(the energy dependence of the tree contribution into the exponent is by now a technical matter, 
as I will discuss below). 

°In fact, the results (H), (^) and (0) have been obtained by making use of the functional technique of ref. 0. I 
will discuss this technique in some detail later. 



Thus, the perturbative calculations of processes with multiparticle final states (both in the 
presence and without instantons) provide important insight into the functional form of the 
cross sections. However, these calculations are reliable only at relatively small n when the cross 
sections are still exponentially suppressed. Attacking the most challeging problem of really 
large n requires the development of novel non-perturbative techniques which, as suggested by 
Eqs. (I) and (|), should be of semiclassical type. Before turning to some of the attempts in this 
direction, let me make a comment concerning the expectations for the behavior of multiparticle 
cross sections, which is based on rather general grounds. 



2.3. Unitarity 

The argument I would like to present is in the spirit of refs. and is based on unitarity 
and ordinary perturbation theory at low energies. It shows that the multi-particle cross sections 
are likely to be small at all energies, except, maybe, for exponentially high ones. Consider 
theory for definiteness (the argument works also for instanton-like processes). The propagator 
in this theory obeys the Kallen-Lehmann relation, which I use at small Q, say Q = 0, 

G{Q = 0)=J ds^^ (8) 

(the argument is easily generalized when subtractions are necessary), where (Ttot{s) is the total 
cross section of the " decay" of one virtual boson with energy ^/s into arbitrary nubmer of real 
bosons. The left hand side of this relation is believed (and in some cases proved) to be a nice 
asymptotic series in A whose finite number, /c -C 1/A, of terms are given by perturbation theory. 
These k terms on the right hand side come from the cross sections of production of k particles 
and less. Since the remaining part of the left hand side is smaller than const ■ A'^^^, the sum of 
contributions of n-particle cross sections with n > k into the right hand side is small, 

g f rfs^ < const ■ A'^+i 

n=k+l •' ^ 

This excludes the possibility that the ri-particle cross sections at n ~ 1/ A are large at energies 
of order m/A. In fact, in view of Eq. (|^), these cross sections are expected to be exponentially 
small, i.e., the exponent F is expected to be negative at any n. 

This argument does not, however, exclude a still very intersting possibility that F tends 
to zero as An — ^ oo and/or e — *• oo. In any case, theoretical understanding of the exponential 
behavior of the cross sections is of considerable importance. 

Equation (||) actually provides a possible way to estimate (or rather put an upper bound 
on) n-particle cross sections 0. The left hand side can be non-perturbatively calculated on a 
lattice and then compared to a few terms in the perturbative expansion of G{Q = 0). The 
difference is then attributed to the contributions of n-particle cross sections at large enough n. 
This approach has already been tested in (1 + l)-dimensional scalar theories 0. 



3. Regular Classical Solutions: from many — > many to few many 

One possible way MM to study the probability of the processes few many is to consider, 
as the first step, the processes many many, i.e.. 



where , the number of incoming particles, is formally of order the inverse coupling 
constant. For finite auin and small a, both the initial and final states contain parametrically 
large number of particles, and it is natural that the scattering process can be described in 
semiclassical terms. The probability of few — > many process may then be obtained in the limit 



One realization of this idea is to consider real classical solutions to Minkowskian field equa- 
tions, i.e., scattering of classical waves. To every classical solution that disperses into free waves 
at t — *• ±cxo one can assign the number of incoming and outgoing particles, both of which are 
naturally of order 1/a. The probability of the scattering of these multiparticle states is not 
suppressed. At given energy one tries to minimize the number of incoming particles under 
the condition that the topological number changes by 1 (instanton-like transitions) or that the 
number of outgoing particles is fixed (processes without instantons). If the minimum number 
of incoming particles tends to zero (in units 1/a) as the total cm. energy approaches some 
Ecr, then few — > many processes are not suppressed exponentially ai E > Ecr (this includes 
more likely possibility that Ecr = C)0, in which case the exponential suppression of few many 
cross sections disappears at asymptotically high energies). In the opposite case when the mini- 
mum number of incoming particles needed to induce the classical transition remains larger than 
const ■ the exponential suppression of few —>■ many persists at all energies, but the actual 
exponent cannot be calculated by studying classical scattering. 

Naturally, obtaining general enough set of analytical classical solutions is possible only in 
a very narrow set of models. In specially designed models @ (false vacuum decay in scalar 
theories with potentials of exponential type in (1 + 1) dimensions) it has been found that the 
number of incoming particles required for the classical instanton-like transition to occur is finite 
in units 1/a at all energies. This provides an example of few —>■ many transitions whose rate is 
exponentially small at arbitrarily high energies. 

In realistic models, this program requires extensive numerical calculations. Most impressive 
results up to now have been obtained in ref. where (3 + l)-dimensional SU(2) Yang-Mills- 
Higgs theory has been studied. It has been found that, indeed, the minimum number of incom- 
ing particles producing classically instanton-like (or rather sphaleron-like) transitions decreases 
with energy. Though the results of ref. are not decisive yet (in these calculations dropped 
by about 30 per cent only), this study suggests that the full scale reahzation of this program 
is quite feasible. A remarkable feature of this study is that it proved "experimentally" that 
distinguishing processes in a trivial vacuum and sphaleron-like transitions is indeed possible, 
and one is really able to study the latter transitions even though topologically trivial processes 
may be much more numerous. 

When the number of incoming particles and total energy are such that the rate is exponen- 



is able, however, to formulate a classical boundary value problem in complex time S To fix 
the number of incoming particles and total energy, one introduces the Legendre conjugate real 
parameters, T and 6. The contour in complex time plane on which the classiacal field equations 
are solved has to start at T = t' + iT/2, where t' = real infty. At this asymptotics the 
boundary condition for (spatial Fourier transform of) the field is 



'm 



0. 





(9) 



where 6k are arbitrary complex functions of k. The contour should end at real t ^ oo, and 
another boundary condition is that the field is real in the future asymptotics. The number of 
the boundary conditions is sufficient for determiming the classical solution, up to translations 



in space-time. The total probability of the instanton-like transitions indeed has the exponential 
form, provided that aE amd ariin are fixed as a — >• 0, 



atot oc exp 



—F{aE, arii. 
a 



The solution to the above boundary value problem determines the exponent, 

-F(aE, aN,n) = ET - nO - 2Sd 
a 

where Sd is the real part of the classical action with Euclidean time convention, and T and 9 
are related to E and rij^ by 

Thus, the study of many — >■ many transitions in the classically forbidden region also reduces to 
the problem of solving classical field equations, but now for complex classical fields in complex 
time plane. Analytical results are again available only in specially designed models S where 
they nicely match those obtained by studying the classical scatteringil. There has been very 
few attempts to find the corresponding classical solutions numerically eIS. The most intersting 
results up to now have been obtained in ref. where false vacuum decay in (p^ theory has been 
studied in (3 + 1) dimensions. The energy range accessible was from Q.^Egph to SEgph] it 
was found that at the highest energy in this interval, the exponential suppression disappears 
at = OAugph, where Egph and ngp^ are energy and number of particles charactristic to 
the sphaleron decay, Esph-,nsph ~ l/a. The results of ref. indicate that the exponential 
suppression does not disappear at high energies at nj„ <^ risph, and this expectation should be 
confirmed or disproved in near future. 

Thus, the two approaches which make use of the idea of many many transitions comple- 
ment each other and are capable to provide a coherent picture of the instanton-like processes at 
high energies. The existing results indicate exponential suppression of few many processes, 
and decisive results are expected to be obtained soon. 



4. Singular Classical Solutions: towards the Generalization of Landau Technique 
to Quantum Field Theory 

It is well known in quantum mechanics of one variable that semiclassical matrix elements can 
be evaluated by making use of singular classical solutions to equations of motion in Euclidean 
time ii. Namely, one considers matrix elements of the type 

{E2\d\Ei) (10) 

where \Ei) and I-E2) are highly excited (semiclassical) states, and O is an operator independent 
of h. The Landau technique states that this matrix element is exponentially small, 

{E2\d\E^) (X e-^l^ 

and the following properties hold: i) the exponential factor is independent of the operator O; 
ii) The exponent S is equal to the truncated action of a singular classical solution in Euclidean 
time, which has energies Ei and E2 before and after the singularity, respectively; this solution 
begins and ends at classical turning points. 



It has been pointed out by several authors Eac2lcllE3'E3 that there exists a remarkable sim- 
ilarity between the matrix elements (0) and amplitudes of multi-particle processes at high 
energies, so the Landau technique may be generalizable to field theory. In fact, this technique 
was successfully used for the evaluation of tree amplitudes at threshold in 0^ theory H and 
also for estimating full amplitudes at threshold in this theory at asymptotically large number 
of final particles e3. 

Indeed, consider, as an example, multi-particle production in theory (no instantons). 
The corresponding amplitudes are 

{n\d\Q) (11) 

where both vacuum and final states may be viewed as semi classical states at n ~ j, and O 
is 0(x) or (j){x)(j){y) (for one-particle or two-particle initial states, respectively); O obviously 
does not depend on A. The similarity between Eqs. (p!0| ) and (|Tl]) is clear. Further support 
to the idea to generalize the Landau technique comes from perturbative calculations which 
show that the exponent for the amplitude near threshold is in fact independent of the operator 
O in Eq. (|11]), at least at exponentiated one loop order for operators like 0(x) and (f){x)(f){y). 
Most notably, singular solutions appear naturally in the calculations of tree cross sections at 
and above threshold. 



4-1- Singular Solutions and Tree Cross Sections 



Let me discuss in more detail the relevance of singular classical solutions to the tree cross sec- 
tions. I am going to show that the tree n-particle cross sections at A ^ with Xn, E/n =fixed, 
are related to singular solutions to Euclidean field equations. The resulting prescription has 
been suggested in ref. ^1 on slightly different grounds. 

To begin with, consider generating function for 1 — >■ n tree amplitudes at n-particle thresh- 
old. It has been shown 11 that it obeys Minkowskian field equation 



and it is homogeneous in space. 



9^0 + + A03 = 
: 4>{t). The relevant solution is 
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(12) 



where Zq is a free parameter. The tree amplitudes at threshold are obtained by writing this 
solution in the following form. 



z{t) 



and expanding in z{t), i.e., 

(^,Pi 
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J 2=0 



In this way one recovers the amplitudes (^. 

Being continued to complex time, the solution (|l^) decays at Im t ^ oo, and has a singu- 
larity at 



— m 

2m 



Xzl 



Viewed as the field configuration in Euclidean time, this solution has flat surface of singularities 
near which 

AT 

where / = Im t — TQ = T — Tois the distance to this surface. The flatness (independence of x) of 
this surface can be viewed as a reflection of the fact that all spatial momenta of outgoing and 
incoming particles are zero. 

This way the singular solutions emerge may appear accidental. To see that this is not so, 
consider 1 —>■ n process at the tree level above ra-particle threshold. Clearly, the amplitudes 
themselves may depend on the correlations between the momenta of the outgoing particles, so 
the form of the amplitude is not expected to be simple. On the other hand, the total (integrated 
over phase space) n-particle cross section depends only on n and total energy E. It is therefore 
natural to look for a semiclassical way of calculating this cross section. 

As the first step, let me consider the amplitude of the production of a coherent state |{6(k)}) 
at given energy E. It is given by the functional integral in the coherent state representation 
with appropriate boundary terms at t = Ty ^ oo, 

({6(k)}|<^(^)|0) = / e^s+Bit,^,) ^14) 

where 0/ = (piTf), 

B = -\j dkb*{k)b*{-k)-^J dkuk(Pf{k)(Pf{-k) + J dkV2^kb*{k)(Pf{-k) (15) 



and 



4){E) = J d^dt 0(x)e-'^*+*^'' 



Here ~ 1/A is the total cm. energy, and P = is the total cm. momentum. The amplitude 
(p^) is determined at the tree level by the saddle point of the exponent, which is a solution to 
the classical equation 

+ m^(i>c + A0|! = 

The boundary conditions for are obtained by varying the fields at t = Ty — >■ +oo and 
t = Ti ^ — oo. One finds that 0c has only positive-frequency part at t — — oo (Feynman 
boundary conditions), 

0c(k,t) oce^'^'=*, t^-oo (16) 
and at t — s> +oo its positive-frequency part is determined by &*(k), 

0c(k, t) = — l^e*^"* + (arbitrary) ■ e'^^'^*, t +oo (17) 



In fact, the energy of the field ([T6| ) is zero, so 0c has only positive- frequency part at t ^ +cxo 
by energy conservation, 

0c(k,t) = ^e-^-* t-^+oo (18) 

The exponent in Eq. ( p!^ turns out to be zero for = 0c, so the tree amplitude has particularly 
simple form, 

{{b{k)mE)\0) = AE{b*) = I d^dtUb*;x)e-''''+'^^ (19) 



The second step is to extract the n-particle cross section from the amphtude (|T^). To do 
this, one notes that, in the spirit of ref. Eil, the generating function for n-particle cross sections 
has the following form, 

S(e; E) = ^J DbDb* exp (| dk 6(k)6*(k)) AEi^b*)AEi^^b) (20) 

where F is the usual flux factor. Indeed, the n-particle amplitude is proportional to the term 
containing a factor 6*(ki) . . . &*(kn) in the expansion of AE{b*) in 6*'s. The expansion of Eq. 
(pop in ^ is thus the expansion in the number of outgoing particles, and integration over b, 
b* performs, as usual, the integration over phase space. Thus one obtains, after a change of 
variables, 

^rZXE) = jj^, DbDb* exp (^-1 J dk 6(k)6*(k)j A^(6*)A^(6) (21) 

where the integration contour in complex ^ plane surrounds the origin. 

I am going to perform the integration in Eq. (^1]) in the saddle point approximation, which 
is valid at n ~ 1/A. However, the integration over b, b* contains zero mode due to space- 
time translations. Also, the expression ( |T9| ) contains exponentially large factor due to large E 
[E ~ 1/ A) . Indeed, let me parametrize 

6*(k) = /5*(k)e'^^*°~*'^'' 

and assume that the set {/5*(k)} obeys four constraints (about which I will have to say more 
later). Then, due to the boundary conditions (p!6D, ([1] 



4>c{b*;x) = (j)c{P*;x + xo) 

so that 

AE{b*) = e'^*«-'P"Mi5(/5*) 
The integration over b, b* becomes the integration over (3, (3* and xq, Xq, i.e.. 



exp 



(Ti^niE) = jj J D(3D(3* dtodt'o d^od^'o 
exp [-n\n^ + iE{to - t'o) - 'iP(xo - x^)] 
/ rfk/3(k)/5*(k)e^^'=(*°-*o)-^''(^°-^o) 



where J contains (5-functions of the constraints and the corresponding Faddeev-Popov deter- 
minant. 

Let me now assume that the constraints are chosen in such a way that AeIP*) does not 
contain exponential factors. Then the integral over [xq + x'q) gives rise to the usual volume 
factor in the cross section (which is cancelled by the flux factor 1/F). The integration over 
remaining variables is of saddle point nature. The saddle point value of (xq — Xq) is zero, so I 
omit this variable in what follows. Denoting 

^ = to - to = -iT (22) 

one finds the effective "action" 

Wtree{T, 6, (3*, (3) = ET - uO - e-^ J dk /3(k)/3*(k)e'^^^ (23) 



which should be extremized with respect to T, 9, (3 and (3* under the constraints to which I 
turn now. 

The constraints to be imposed on /3*(k) should fix the translational invariance in space and 
time. Let me impose them implicitly by requiring that Ae{(3*) does not contain exponential 
factors in spite of the fact that E is of order 1/A. One expects that relevant classical solutions 
are generalizations of Eq. (p!2D . Indeed, due to Eqs. ( |T6|) and ( |T8|) the soultions should decay 




as Im t — i> +00. In other words, the solutions considered in Euclidean space-time should 
vanish as r +oo, where r is the Euclidean time coordinate. Since there are no instantons, 
the solutions must be singular in Euclidean space-time. Generically, the points of singularities 
form a three-dimensional surface r = r(x), near which the field behaves like 

(24) 

where / is the distance to this surface (cf. Eq. (0)). The right hand side of Eq. (|19D for such a 
solution behaves like exp(ii^rm + iP^m) where and Xm are the Euclidean coordinates of the 
singular point closest to the Minkowski time axis (r^ < for the solution be smooth everywhere 
in Minkowski space-time). For the right hand side of Eq. (|19D to contain no exponential factors, 
one requires 0, = 0. In other words, the singularity surface in Euclidean space-time 

should touch from below the surface r = at the origin. This constraint fixes the translational 
invariance (both in Minkowskian and Euclidean space-time). 

So, the exponent of the tree cross section is obtained in the following way. One considers 
Euclidean classical solutions 0c(x, r) which decay at r — > +oo and have singularities on surfaces 
T = Ts{x.). For any given singularity surface there exists at most discrete set of such solutions: 
indeed, one may first consider a well defined boundary value problem [ (p ^ at t —>■ +oo; 
(f){T = Ts(x)) = A] and then send A — > oo. The behavior of the solution at large r, 

0(k, r) = ^e-- (25) 

determines the functions /3*(k) for a given surface r^. Then one finds an extremum of the right 
hand side of Eq. (^31) among all surfaces obeying the constraint 

r,(x = 0) = 

and having rs(x 7^ 0) < 0. Finally, the extremum with respect to T and 6 of the resulting 
W{T, 9) determines T and 9 in terms of E and ra, and the cross section is 

oc exp {Wtree[T{E, u), 9{E, n)]) 

This prescription coincides with one advocated in ref. @. 

The dependence on A of the tree ra-particle cross section is, in fact, obvious from Feynman 
graphs, cr^!^^ oc A". This dependence is easily restored within the above prescription by the 
following scaling argument. Since the integral 

J rfk /5(k)/3*(k)e'^'=^ (26) 

as well as the constraints on the surface of singularities, do not depend on 9, the extremum 
value of this integral depends only on T, 



/ 



dk /?(k)/5*(k)e'^^-^ = \n{T) 
A 



(the dependence on A comes about from even more simple scaling), where N{T) should be 
determined by solving the classical problem oulined above. Therefore, the functional form of 

Wtree IS 

WtreeiT, 6) = ET - uO - e-'jN{T) 

Then the extremum with respect to 6' is at ^ = — In and the value of Wtree at this extremum 
is 

Wtree{T) = ET + u ln{Xn) - n In N{T) (27) 
After extremizing this expression with respect to T, one finds the functional form of Wtree{E, n), 

WtreeiE, u) = u In(An) + n^(e) 

where, as before, e = E/n is the average energy per an outgoing particle. Therefore, the 
functional form of the tree cross section is 

aiZl ^ oc n!A"e"*(^) (28) 

which has correct dependence on A. Note that the tree cross section increases with n at given 
e, i.e., it indeed hits the unitarity bound at large enough n. 

Let me also point out that the extremum of the integral ( p6D is naturally expected to be 
a minimum; the perturbative calculations suggest that the extremum of W{T), Eq. (pTj), is 
also a minimum. [The latter observation is in accord with the fact that the integration over 
the original variable (to — ^q) runs along the Minkowski time axis, and the extremum of W is 
a maximum along this line. Switching to T according to Eq. (|2^) converts the extremum into 
a minimum.] Therefore, it is possible to formulate a variational procedure for obtaining the 
lower bound on \l/(e) entering Eq. (^); this procedure has been applied to 0^ theory in ref. 0. 

Thus, there exists a well defined procedure that relates the exponent for the tree n-particle 
cross section to singular solutions of classical field equations in Euclidean space-time. The 
existing perturbative results have been reproduced within this approach @, and the evaluation 
of the behavior of tree cross sections at all energies becomes a matter of numerical calculations. 

Loops 

As discussed in sect. 2.1, loop effects become important at An ~ 1, and the correct behavior 
of the n-particle cross section is determined by diagrams with arbitrary number of loops. One 
may argue, however, that the complete exponent for the cross section, including loop effects, 
is again related to singular classical solutions of the field equations. One argument is close in 
spirit to the approach outlined in sect. 2. 2. Consider, instead of the matrix element (14), the 
following matrix element, 

A{j) = {{bik)}\PEexp(^J rfxj(a;)0(x)) |0) (29) 

where Pe is the projector onto the subspace of total energy E, and 

j{x) =jo6{x) 

is an external source. The projection by Pe may be written in the integral form 0, so that 

A{j) = J dr] e*^''({6(k)e'"*''}|e^»'^(°)|0) 



and the matrix element (^) has the functional integral representation similar to Eq. (|T^) . 
It is then straightforward to write the double functional integral representation for the "cross 



section" of the creation of n particles by the external source j; in particular, one makes use of 
the trick (^Of ). The form of this double functional integral is, roughly speaking. 



aj^n = J D(l)D(f)' DbDb* dT]d^ 



exp - iS[(l)'] + jo0(O) + jo0'(O) +iEr] + n\n^ + ...) (30) 

where dots stand for boundary terms like ( |15|) and terms proportional to bb*; all these terms 
are bilinear in [b, b*, (j), (j)'). It is clear that at 

1 

the integral (0) is of saddle point nature, and its value is determined by a solution to the 
classical field equations in the presence of the source, which obeys the boundary conditions 



similar to Eqs. (|T^) and (|T^. The "cross section" has then the following functional form, 

(31) 



(Tj^n exp 



nF ( An, — , Ajo 

n 



The amplitudes few — n are obtained from the matrix element (^) when jo is small. Therefore, 
one argues that the exponent of the full cross section few n is obtained by taking the limit 
■^Jo ^ in F {\n, E/n, XJq). In this limit the saddle point configuration obeys sourceless 
classical field equation. Clearly, this configuration must be singular, otherwise the energy of 
this configuration would be the same at t — >■ —oo and t —>■ +oo, while these energies should be 
equal to and E, respectively. The nature of the singularity is to be determined by the above 
limiting process. 

This idea has been elaborated in ref. where the classical boundary value problem for 
fields with subsequent extremisation over their singularity surfaces has been formulated. In the 
limit of small An, when the tree cross section is a good approximation, this procedure reduces 
to one discussed above. The exponentiated one loop correction has been also reproduced by 
this technique So, known results summarized in Eq. indeed follow from the semiclassical 
treatment that generalizes the Landau technique to field theory. It remains to be understood 
what kind of singularity of the calssical field is relevant in the general case, and whether this 
technique is suitable for numerical studies. 



5. Summary 

To summarize, few many cross sections cannot, in general, be understood within per- 
turbation theory, irrespectively of whether one considers instanton-like processes or processes 
in the trivial vacuum sector. The functional form of the cross sections in both cases strongly 
suggests the possibility of semiclassical treatment of these processes. Several semiclassical ap- 
proaches have been recently developed and tested analytically and numerically. The decisive 
results have not been obtained yet, but they are expected to appear in near future. 

This work is supported in part by INTAS grant 94-2352 and International Science Founda- 
tion grant MKT300. 
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